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Abstract. Addition of higher nonlinear terms to the well known integrable nonlinear 
Schrodinger (NLS) equations, keeping the same linear dispersion (LD) usually makes the 
jjy ■ system nonintegrable. We present a systematic method through a novel Eckhaus-Kundu 

hierarchy, which can generate higher nonlinearities in the NLS and derivative NLS equations 
preserving their integrability. Moreover, similar nonlinear integrable extensions can be made 
again in a hierarchical way for each of the equations in the known integrable NLS and 
derivative NLS hierarchies with higher order LD, without changing their LD. 
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^ ; 1 Introduction 



Nonlinear Schrodinger equation (NLSE): 



Ft ' iq t + q xx - 2a\q\ 2 q = 

> 

as well as derivative NLSE (DNLSE) 



iqt + Qxx ~ ia(\q\ 2 q) x = (2) 

arising in different physical context [1] , represent well known integrable systems. The inte- 
grability of such systems as well as the stability of soliton solutions inherent to such equations 
are believed to be due to a fine balance between their linear dispersive and nonlinear collapsing 
terms. Therefore, in some physical situations, which demand addition of higher nonlinear terms 
to (1) or (2) [2, 3, 4], this balance apparently is lost and the system turns into a nonintegrable 
one, not allowing analytic solutions [5]. On the other hand, together with the well known NLSE 
and DNLSE (1) - (2) there exists a tower of equations in their integrable hierarchies correspond- 
ing to higher conserved charges, where increasingly higher nonlinear terms do arise. However at 
the same time higher order linear dispersive terms also appear in these integrable equations, ap- 
parently to compensate for the higher nonlinearities and for restoring the balance. For example, 
the next equation in the NLSE hierarchy (evolving with time £3) has the form 

qt-s + qxxx + Qcr 2 \q\ 2 q x = (3) 

with higher nonlinearity ((/l 2 ^, which however has a compensating higher order linear dispersive 
term: q xxx . 
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Nevertheless, in physical systems quite often one needs to add higher nonlinear terms (includ- 
ing nonlinear derivative terms) without increasing the linear dispersion (LD) of the system. An 
excellent example of such a system is the Johnson's equation [2], which was derived to solve an 
important hydrodynamic problem for analyzing Stokes instability in fluid flow near the critical 
value of kh. This equation which is uniformly valid for any kh is given in the form: 

A t - axA xx - a 2 \A\ 2 A + a 3 |A| 4 ,4 + ia 4 \A\ 2 A x - ia 5 (\A\ 2 ) x A - a 6 9 t A = 0, (4) 

where ai are real numbers and 6 X = <5|^4| 2 . Evidently this physical system contains the same 
linear dispersive (LD) term a\A xx as in the NLSE and DNLSE together with their associated 
nonlinear terms, whereas it also has additional 5-th order nonlinearity along with nonlinear 
terms containing derivatives (described by the terms with coefficients 03,04 as well as a§). As 
shown in [2], at criticality: kh = 1.363 and all ai have known values with 0,2 = 0. 

Another example of such a physical system is the Benney's hydrodynamic equation [3] de- 
scribing long and short wave interaction. Therefore it is natural to ask whether one can add 
higher nonlinearity to the original integrable NLS and DNLS equations (l)-(2) without changing 
their LD, and at the same time preserving their integrability. Moreover, we may enquire whether 
it is possible to do this in a systematic and hierarchical way. 

We focus on these intriguing questions and find that there exists an integrable hierarchy of 
Eckhaus-Kundu equations [6, 7], extending the NLS and DNLS equations to higher nonlinearities 
and thus it provides a conclusive answer to the above posed questions. Moreover, such a hierarchy 
with the addition of nonlinear terms in a recursive way without changing their LD and the 
integrability can also be constructed for each of the equations with higher order LD in the 
known integrable NLSE and DNLSE hierarchies. 

2 Higher nonlinear NLSE and DNLSE 

It has been shown [6] that, under a nonlinear transformation of the field q — > Q = qe~ lS with 
arbitrary gauge function 6 the NLSE yields an integrable higher nonlinear equation 

iQt + Qxx - 2a\Q\ 2 Q - (0 t + 9 2 x - i8 xx )Q + 2iO x Q x = (5) 

The DNLSE is extended similarly under the same field transformation to 

iQt + Qxx - ia(\Q\ 2 Q) x - {0 t + 9l- i6 xx )Q + 9 x (2iQ x + a\Q\ 2 Q) = (6) 

It is interesting to note that the generation of such new integrable equations can be linked to 
the gauge transformation of the corresponding Lax equations : & x = 3>t = of the 

NLSE and DNLSE to the new systems $ x = U$, $ t = V$, with the gauge transformed Lax 
operators U = hUh~ l + h x h _1 and V = hVh^ 1 + hfh -1 , where the gauge matrix h G U(l) 
is given by h = e t6r7s , with arbitrary gauge field 9(x,t). Note that though both the equations 
(5), (6) have been studied quite extensively [8, 9, 10], the investigations were confined mostly to 
the particular functional choice for 9 as 9 = 5 f x \Q(x')\ 2 dx', with a real parameter 5. Under 
this choice the NLSE (1) is extended to the Eckhaus-Kundu (EK) equation [6, 7, 11, 12] 

iQt + Qxx - 2a\Q\ 2 Q + 5 2 \Q\ 4 Q + 2i5{\Q\ 2 ) x Q = (7) 

while the DNLSE (2) turns into a similar equation 

iQt + Qxx - ia(\Q\ 2 Q) x + 5 -(25 - a)\Q\ 4 Q + 2iS{\Q\ 2 ) x Q = (8) 

Note that in both the above equations higher nonlinear terms including derivative and 5-th power 
nonlinearity are added to the NLSE and DNLSE, without changing the linear dispersive term 
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and without spoiling the integrability of the system. It was detected [6] that, for the parameter 
choice a = 28 and a = 8 , (8) turns into two other well known equations, e.g. Chen-Lie-Liu [13] 
and Gerdjikov-Ivanov equations [14]. 

We intend to generalize now this concept and find the general form for additional nonlinear 
terms, for both NLSE and DNLSE, such that they can be included without altering the linear 
dispersion and the integrability. We thus find a novel hierarchy of EK equations which provides 
a systematic way for such construction. We observe first that in equation (5) the gauge function 
B should be chosen through its x and t derivatives as B x = 8\q\ 2 , Bt = iS(q*q x — qq x ), f° r obtaining 
higher nonlinear NLSE (7), while one should choose in (6), B x = 8\q\ 2 , Bt = i8(q*q x — qq x ) + 
|«(5|g| 4 , to get the extended DNLSE (8). It is important to note that, in both these cases the 
necessary condition (B x )t = (0i) x must hold, which follows here from the validity of the NLSE 
and DNLSE themselves for the field function q(t, x). 

Our idea, therefore, is to widen this choice of B for finding a hierarchical rule, by using the 
conservation relation ip t + I x = 0. Such integrable systems, as is well known, possess infinite 
number of conserved quantities in involution, and we can make a systematic and recursive choice 
for B x n ^ = p( n "> and B^ = il^ n \ n = 0, 1, 2, . . . for constructing our hierarchy. 

Interestingly, we can apply the idea of this construction to each of the higher equations in the 
NLSE and DNLSE hierarchies with higher order LD to generate a new integrable EK hierarchy 
with higher nonlinear terms, but with the same LD. 

To remove any possible confusion, we should remark here that the integrable EK hierarchies 
we discover here for the NLSE and DNLSE, together with the known equations (7, 8) as well 
as the well known Chen-Lie-Liu and Gerdjikov-Ivanov equations etc. do not appear in the 
celebrated classification list of Mikhailov et al [15], since it excludes equations which are related 
to the enlisted equations through gauge or invertible nonlinear transformations. 



3 Recursive formulas and general framework for conservation 
law 

We first present a general framework for deriving the required conservation rule for both NLSE 
and DNLSE. We start with the linear set of Lax equations: $ x = ® t = V<& with $ = 
($i,$2) and eliminate $2 from both these equations. That would yield from the first Lax 
equation a second order linear equation for <J>i: 

*ixx = (^P)*!* + (U llx + U u + U 12 U 21 - ^pC/n)$i, (9) 

where Uij are matrix elements of the space-Lax operator U. Representing now $1 = e u ^ x +<P 
the linear equation (9) turns into a Riccati equation for (p x = v(X): 

U l2 {-^-) x + 2U ll v + v 2 = U 12 U 2U (10) 

which by expanding v{X) in powers of A -1 generates an infinite set of conserved densities c n , n = 
0, 1, 2, — On the other hand from the time-Lax equation we derive 4>t = Vu + V\ 2 while the 

space-Lax equation gives ^ = yielding <fi t = Vu + j^<f) x . After taking x-derivative this 
gives the crucial relation 

vt = (Vu + ^ v) x (11) 

Note that we can derive now the conservation law and hence the hierarchical expressions for the 
gauge field B from (11) for concrete systems of NLSE and DNLSE by inserting the corresponding 
expressions for the Lax matrix elements Uij, Vij. 
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4 Higher nonlinear hierarchy from NLSE 

We concentrate first on the NLSE system for deriving its integrable hierarchy of higher nonlinear 
equations with the same linear dispersive term. For NLSE the conservation relations were 
established explicitly in [16]. Using the known form for the Lax matrices of the NLSE: 

Uu = -iX, U12 = \foq , U 2 i = \/aq* 
and V u = -i(2X 2 + a\q\ 2 ), V 12 = V^(2Xq + iq x ), (12) 

we derive the corresponding Riccati equation from (10) as 

2i\ v = u 2 - a\q\ 2 + q(^) x . (13) 

Since v vanishes at |A| — > oo, expanding it in spectral parameter v = J2 n =o (2ixfcn ■> we ^ n( ^ a 
recursion relation from (13) for the densities of the conserved quantities: 

n-l 

Cn+i = q{—)x + ^2 c fcCn-fc-i, for n > 1 (14) 
9 k=0 

with Co = — cr\q\ 2 , c\ = —aqq*, c 2 = —a{qq* xx — o"|<?| 4 ), etc. This gives a methodical way for 
evaluating the infinite set of conserved quantities C n = J c n dx,n = 0,1,2,.... Since these 
conserved quantities are related to the NLSE, which is an integrable system in the Liouville 
sense, all of them must be in involution : {C n ,C m }pB = [16, 17]. This may be checked 
by using the fundamental PB relation {q(x), q*(y)}pB = 8(x — y) and the boundary condition 
lim x ^± 00 \q\ — > 0. We focus now on crucial relation (11), which using the relevant Lax matrix 
elements for the NLSE derives the conservation relation 

v t = -i(a\q\ 2 + (2i\-^)v) x (15) 

Expanding further u(X) in parameter A, we get the hierarchy of relations for the densities of 
conserved quantities 

C n t = i ( -Cn+l + — Cn J , 71 = 0,1,2,... (16) 

where all densities c n can be evaluated from the recurrence relation (14). Choosing therefore 
p*™) = —\c n and /( n ) = ^(c n+ i — yc ra ) we can derive finally the required conservation law 

=0,n = 0,l,2,... (17) 
Considering few starting values n = 0, 1, . . . we can easily evaluate their explicit forms as 

p (0) = | g |2, J(0) = _ {qq * _ p (l) = qq ^ Jd) = ( | fe |2 _ qqlx + ^1^4^ (lg) 

etc. Therefore we conclude that we can have a series of choices for the gauge function 9 given 
by 

6^=5[il^\ = 5 -[ l { Cn+1 - q J- Cn )] = ^(c n+1 -^c n + l 5 -c 2 n % (19) 
a q a Q a 

for each values of n = 0, 1, 2, . . ., where [x] = x + x* (or i(x — x*)) indicates nontrivial real valued 
combination of x together with convenient normalization by some constants. For example, one 
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should evaluate the term appearing in the above expression as [i q Cn] = i(^c n — c*). Note 
that by inserting this choice for 9^ in higher order equation (5) one can generate a novel 
hierarchy of EK equations for n = 0, 1,2, . . .. At the bottom of this hierarchy, as we see from 
(18), lies the simplest equation (7). However it is important to note, that since the gauge field 9 
in our construction should be a real function, when it is chosen through conserved densities c n , 
which in general can be complex valued as evident from (14), real combination of its expression 
must be taken with proper care, along with any convenient normalization. Such a manipulation 
with c n obviously does not affect the conservation law (16), since it is a linear equation in c n . 
For example, for n = when Co is real we can take 9 X °^ = 5\q\ 2 , however for n = 1, when c\ is 
complex valued, we must have 9 x 1 \q) = — f [ci] = —iS(q*q x — qq*) to make it real. Similarly we 

get for n = 2, 9 x 2 \q) = -£[c 2 ] = 5(q*q xx + qq* xx - 2a\q\ 4 ). 

The simplest and the lowest in the EK hierarchy with n = is given explicitly by (7). 
We therefore present the details for the next higher nonlinear equation with n = 1 in the EK 
hierarchy using (18), while the other equations with n > 2 can be derived in a similar way. 
Following the above formulation we get for n = 1: 

eg\q)=i8T(q), 9^(q) = -5{E(q)+2a\q\ A ) (20) 

where 



T(q) = (q*q x -qq* x ), E(q) = 2q* x q x - (qq* x + q*q xx ). (21) 

The consistency of (20) can also be checked independently by using the NLSE for the field 
q. However, since the equation for the transformed field q — > Q with q = Qe i9 ^ is of 
our interest, we have to rewrite all the above expressions completely in terms of the new 
field, by expressing the functions 9 x (q),9t(q) through Q. For this purpose we derive q*q x = 
(Q* e -iO(Q)j ((q x + i9 x (Q)Q)e ie ^) = (Q*Q X ) + i9 x {Q)\Q\ 2 , which evaluates 

T(q) = T(Q) + 2i9 x l \Q)\Q\\ E(q) = E{Q)+A(-i9 x l \Q)T{Q) + 9 x l \Qf\Q\ 2 ) (22) 
where T(Q) = (Q*Q X - QQ* X ), E(Q) = 2\Q X \ 2 - (Q*Q XX + QQ* XX ) (23) 

resulting 9 ( x \q) = i5(Q*Q x - QQ* X ) - 259 { x 1] {Q)\Q\ 2 which gives 

^ )(Q) = i M(Q) T(g) ' M ^ 1 + 2( ^| 2 ( 24 ) 
and similarly 

9i 1) (Q) = 5(E(Q) + 2a\Q\ 4 + 4(-ie£\Q)T(Q) + ^ 1} (Q) 2 |Q| 2 )) (25) 

where T(Q),E(Q) are as defined in (23). 

As a result we derive the next higher nonlinear equation with n = 1 in the EK hierarchy from 
the NLSE as 

iQt + Qxx~2a\Q\ 2 Q -5((E(Q)+2a\Q\ 4 + 1 ^T x (Q))Q + 2 J l^T(Q)Q x ) 

with E(Q), M(Q),T(Q) as defined above. Note that this is an integrable equation with higher 
nonlinearities, but with the same second order LD as in the NLSE. Notice that (26) contains 
an additional coupling constant 5 entering in different powers. Therefore we may simplify this 
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equation by considering S to be small and ignoring all terms with higher powers in 6, which 
yields 

iQt + Qx X - 2a\Q\ 2 Q + 25(-a\Q\ 4 Q - Q*(Q X ) 2 + Q* XX Q 2 ) = 0. (27) 

However it should be noted that though the simpler equation (27) as such is a nonintegrable 
system, it might be meaningful for physical applications and useful information can be extracted 
for it from its integrable variant (26) through limiting procedure. 

The general form of the EK hierarchy can be generated similarly in the same form (5) where 
the series of gauge functions 9[ n \o x n \n = 0, 1, 2, . . . are chosen as (19). This would construct 
finally a new hierarchy of integrable EK equations emerging from the NLSE in the form 

iQt + Q xx ~2a\Q\ 2 Q +-(([i(-c n+1 + ^-c n )]-i[c nx })Q-2i[c n }Q x ) = 0, (28) 

where the real-value sign [•] should be properly evaluated in explicit calculations. The conserved 
densities c n = c(Q) n ,n = 0, 1,2, appearing in equations (28) are to be obtained systematically 
from (14) and then expressed consistently in terms of the field Q, as we have demonstrated in the 
case n = 1. We stress again that the integrable hierarchy of EK equations (28), contains higher 
and higher nonlinearities (including nonlinear dispersive terms), but has the same second order 
LD as in the original NLSE. Note that the density of the conserved quantities, as evident from 
(14), introduces an additional derivative resulting the appearance of higher nonlinear dispersive 
terms at each higher step. However this should not be confused with the linear dispersive term 
Q xx , which remains the same for the whole EK hierarchy, which in fact is our aim. 



4.1 LD preserving EK hierarchy from NLSE hierarchy 

Recall that the well known integrable NLSE hierarchy is generated by the Lax pair (U, V^)nlse, 
k = 1, 2, 3, ... by the flatness condition Ut k — V x k ^ + [U, V^] = 0, with the kth equation having 
the kth order LD . k = 2 corresponds to the NLSE (1), while k = 3 yields the integrable equa- 
tion (3). It should be noted also that while the space-Lax operator U is the same as that of 
the NLSE, the time-Lax operator V^ k \ k = 1, 2, 3, . . . is different for different equations in this 
hierarchy with higher time tk- Nevertheless, it is interesting that, the idea we have developed 
above for generating nonlinear EK hierarchy from the NLSE preserving its LD, can be success- 
fully implemented to each of the k equations in the known NLSE hierarchy. Moreover, since 
the Riccati equation (13) and consequently the recurrence relation (14) yielding the conserved 
quantities depend only on the Lax operator U, they should be the same for all equations in 
the hierarchy with any k. Therefore one can make the same choice (19) for the gauge function 

x n \n = 0, 1, However since the time evolution (15) with respect to higher time tk would 

(k) (k) 

depend explicitly on the matrix elements , V 12 , k = 1,2,3,... as evident from (11), this 

(n) 

would be different for different k, leading to new choices for the 6\ , obeying the consistency 

(9 x n) ) tk = (0t n) ) x . Thus the kth equat ion in the NLSE hierarchy (with kth order LD) would gen- 
erate EK type nonlinear hierarchy with n = 0, 1, . . . , preserving its linear (kth order) dispersive 
term and the integrability. 

Since we have presented above in detail the k = 2 case, i.e the NLSE, we give here in brief 
the construction of the EK hierarchy for the next equation (3) with 3rd-order LD. It is easy to 
check that (3) can be extended through gauge transformation to the hierarchical form: 

Qt 3 + Qxxx + 6a 2 \Q\ 2 Q x + 3ie x n ^Q xx + 3(ie x n J -(9^) 2 )Q X 

- (-^ + ™ { M n J - + ^ f + Gi*% n) \Q\ 2 )Q = 0. (29) 
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For deriving now v(X)t 3 and consequently 0j£' we have to use , as mentioned above, the relevant 
elements of the corresponding time-Lax operator given by [17]: 

4 3 ) = 4iX 3 + (2iX\q\ 2 -(qq* x -q*q x ))a, 

V x f = (-4X 2 q + 2iXq x -2a\q\ 2 q)^. (30) 

Inserting consistent 6^\9 x n \ in (29) we can generate another new EK type integrable higher 
nonlinear hierarchy for n = 0, 1, . . . with the same linear dispersive term Q xxx . 

We however will not derive here the general form of this hierarchy; instead present the simplest 
and the lowest equation with n = 0. We find that the corresponding gauge function in this 
case can be chosen as 

4°) = 5\Q\ 2 , 9^ = 5 (Q*Q XX + QQ* X - Q* X Q X - 3a 2 \Q\ 4 - 35 2 \Q\ 6 + 3iS\Q\ 2 (Q*Q x - QQ* X )) 

(31) 

Inserting (31) in hierarchy (29) we reduce it to its lowest order equation 

Qt 3 + Qxxx + Ga 2 \Q\ 2 Q x + 6{3i\Q\ 2 Q xx + 3(i(\Q\ 2 ) x - 5\Q\ 4 )Q X + (i(\Q\ 2 ) xx - \s(\Q\%)Q 

+ i(Q*Qxx + QQ* XX - QZQ X + 3a 2 \Q\ 4 - 45 2 \Q\ 6 + 3i5\Q\ 2 (Q*Q x - QQ* X ))Q} = 0. (32) 

Note that this is an integrable equation with higher nonlinearities up to 7-th power in the field 
as well as with nonlinear dispersive terms, but having the same LD term Q xxx as in (3). 



5 Higher nonlinear hierarchy from DNLSE 

The steps formulated above may be applied now for finding the nonlinear integrable hierarchy 
from the DNLSE with the same 2nd-order LD as in the original equation. For this purpose 
we shall use the general form of the conservation relation presented in (10)-(11), but customize 
them for this particular case by considering Lax operators [18] associated with the DNLSE: 

■ A2 3 * > i-( q* 
U = -t—a 3 + -X^[ y n 
4 2 \ q 

(y. f * \ 1 

V = i(— - ^X 2 \q\ 2 )a 3 + i ( g 9 Q J , g = -\^(-\ 2 q + 2iq x + 2a\q\ 2 q) (33) 

Note however that unlike the NLSE the required conservation relations for the DNLSE are not 
readily available in the literature in explicit form and have to be derived carefully by inserting 
relevant expressions from (33). This gives from (9) the linear equation for $i as 

*i« = + (-y|j " M^M 2 + +i£|f)*i, M - A 2 (34) 

yielding in turn the Riccati equation for the DNLS system: 

9*(^)z + ^ 2 -^ + /^M 2 = (35) 

Using now the consistent expansion in the spectral parameter: v = iJ2 n =offi> we § e1] the 
recursion relation 

• n j 

-c n+1 = q*(-^) x +i^2c k c n - k , for n > 0, c = --a\q\ 2 , (36) 
q k=0 
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which systematically generates all higher conserved densities starting from 

Cl = a(iq*q x + \a\q\ A ), c 2 = 2a(q*q xx - 2ia\q\ 2 (q*q x + -qq*) - ^-\qf), etc. 



The infinite set of conserved quantities C n = J c n dx,n = 0, 1, . . . for the DNLSE, which is a 
completely integrable system, must be in involution. This can be checked directly by using its 
fundamental PB structure: {q(x), q*{y)}pB = S x (x — y) and the vanishing boundary condition 
for the field: \q\ —>■ 0, at x —>■ ±oo. Using further the Lax matrix elements (33) we derive from 
(11) the crucial conservation law 



v t = ^-|^ + (-|-i^ + a | g |>j . (37) 
Expanding through spectral parameter [i yields the important relation 

1 i / -1x , „,|„|2\„ A „ n i „ - l„.l„|2 



c nt = y--c n+l + (-i-^ + a\q\ )c n j , n = 0, 1, . . . , c = --a\q\. (38) 

For deriving therefore the higher nonlinear hierarchy from the DNLSE we can use the general 
form (6) with the choice for the gauge function 



{ x n) = % (n) ] 



6, 



c 



a 



■n\ i 



l n) = ^ n) ] = -[-k +1 + H4 + a| (? | 2 ) Cn ] 
(y. — Q 

= -hk + i + (-if +a |Q|2 )Cn _f c 2 ]) n = ,l,2,..., (39) 
(y. _■ 

where the real-value sign [■] has the same meaning as explained above. For example, for n = 
with real cq we have 6 X °^ = 5\q\ 2 , while for n = 1 with complex valued c\ we must take 9^ = 
$(i(q*Qx — QQ X ) + ct| , etc. with the similar reason holding also for the choice of 6[ n \ 

This derives finally the new hierarchy of integrable higher nonlinear equations from the 
DNLSE as 



't 



iQt + Qxx-ia(\QrQ) 



X 



~ i (([-^ C ™ +1 " ^ C «] " ilCnxVQ ~ 2i \Cn]Q^J = 0, (40) 

where again all real-valued expressions should be evaluated, carefully and c n — c(Q^)n 

are to be 

obtained recursively from (36), expressed explicitly through transformed field Q. Thus the tower 
of equations (40) with n = 0, 1, 2, . . . would represent a novel EK hierarchy, at the bottom of 
which with n = giving 

= 5\q\\ 6^ = 5{t{q*q x - qq* x ) + |«|g| 4 ), 

lies equation (8). 

It is worth noting that in analogy with the known reductions [6], The parameter choice 
a = 26 in (40) would yield a new Chen-Lie-Liu type hierarchy, while a = 5 would generate 
another Gerdjikov-Ivanov type hierarchy. All these equations however have the same LD given 
by Q xx 

For demonstrating the novelty of hierarchy (40), we take up n = 1, giving the next new 
equation with 

e£\q) = 5(iT(q)+a\q\ 4 ), 6^ (q) = 5{E(q) + 3iaT(q)\q\ 2 + 2a 2 \qf), (41) 
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where 

T(q) = q*q x - qq* x , E{q) = 2\q x \ 2 - (q*q xx + qq* xx ). 

For further application, as was performed also in the NLSE case, we have to express the gauge 
fields o£\q), 9^\q) given by (41) in terms of the transformed field Q and rewrite 

T(q) = T(Q) + 2^)|Q| 2 and E(q) = E(Q) + 4(-i6^ (Q)T(Q) + ^)(Q) 2 |Q| 2 ). 

Note that though the expressions for 9 x l \q) and ( t 1 \q) differ for the cases of NLSE and DNLSE, 
the entries T(q) and E(q) appearing in them have the same form for both these cases. Using 
the above result after some algebra we arrive at the required expressions 



" (Q) = Mi ra+a|Q|4)i 

eP{Q) = S(E(Q) + 3iaT(Q)\Q\ 2 ~ ^{Q)^ + 2a 2 \Q\ 6 + 4(-^W(Q)r(Q) + 6 X 1] {Q) 2 \Q\ 2 )) 

(42) 

where 

M(Q) = 1 - 25\Q\ 2 , T{Q) = Q*Q X -QQ* X and E(Q) = 2\Q X \ 2 - {Q*Q XX + QQ* XX ) (43) 

By inserting (42) with (43) in the general form (6) we derive finally for n = 1 the next higher 
nonlinear equation in the EK hierarchy linked to the DNLSE as 

iQt + ^-ia(|Q| 2 Q) a; -(^ {1) + (^ 1) ) 2 -4l ) )Q + ^ 1) (2^ + a|Q| 2 Q) = (44) 

where expressions for O^iQ) and 6 ( x ] (Q) expressed through T(Q),E(Q) as in (42-43) together 
with 6 x 1 x \q) by extracting the x-derivative are to be inserted in equation (44) to get its explicit 
higher nonlinear form. As we see this integrable equation contains additional higher nonlinear 
terms upto 7th power nonlinearity together with many nonlinear dispersive terms, though with 
the same second-order LD: Q xx as in the original DNLSE. . 

5.1 LD preserving EK hierarchy from DNLSE hierarchy 

The integrable hierarchy of the DNLSE is generated similar to the NLSE by the associated Lax 

pair (U, V^) dnlse, k = 1, 2, 3, The space-Lax operator U corresponds to the DNLSE (33), 

while the time-Lax operator V^ k \ k = 1, 2, 3, . . . is different for different higher time tk, yielding 
through the flatness condition the k-th equation in this hierarchy with k-th order LD. k = 2 
corresponds to the DNLSE (2), while k = 3 yields the next integrable equation [18] 

3 

qt :i ~ q xxx + 3ia(\q\ 2 q x ) x + -a 2 (\q\ 4 q) x = (45) 

The above scheme we have implemented to the DNLSE (with k = 2) can also be applied to each 
of the k > 2 equations for generating integrable EK type hierarchy preserving its kth order LD. 
Moreover, the Riccati equation (35), the determining relation (36) for the conserved quantities 
and consequently the choice for the gauge function 6 x n \n = 0, 1, . . . as given in (39), would be 
the same for all equations in this hierarchy, since they all depend on the space-Lax operator U 
only. However the relations like (37, 38) and consequently the consistent expression for 
which are linked to the time evolution should be calculated using the matrix elements of the 
time-Lax operator V^ k \ k = 1, 2, 3, 
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Thus we can construct in principle higher nonlinear hierarchy with n = 0, 1, . . ., from any of 
the kth equation in the DNLSE hierarchy, preserving its kth order LD and the integrability. In 
practice however each case with higher k would be more complicated due to more complicated 
structure of its time-Lax operator V^ k \ Since we have presented in detail the DNLSE case given 
by k = 2, we report here in brief the construction from the next equation (45) obtained for 
k = 3, which yields 

Qt, - Q xxx + 3ia(\Q\ 2 Q x ) x + l<r 2 ((\Q\ 4 Q) x + iei n WQ) 

- 3i9 x n) Q xx - 3(i0 x n J - {e x n) f)Q x - 3a((\Q\ 2 9^Q) x + 9 x n) \Q\ 2 (Q x ~ ^Q)) 

+ (idff + se^9 x n J-ie x t + i(e^f)Q = o. (46) 

By inserting 9 x n ^ from (39) and deriving the consistent 9^ we can generate from (46) the EK 
type integrable hierarchy with n = 0, 1, . . . , preserving its 3rd-order LD. 

The simplest and the lowest equation of this hierarchy with n = may be given explicitly by 
the choice of the gauge function: 

X 0) = 6\Q\ 2 , 

0<M = 5 (q*Q xx + QQ* XX - Q* X Q X + M(a + 5)\Q\ 2 (Q*Q X - QQ* X ) - 3(25a + 5 2 - \a 2 )\Qf 

(47) 

consistency of which can be checked by using the originating equation (45). Inserting (47) in 
(46) one would clearly get an integrable equation having nonlinearity up to 7-th order in the 
field including its derivatives, but with the linear dispersive term Q xxx as in (45). 

6 Concluding remarks 

Our construction of new EK type hierarchies of integrable equations with higher and higher 
nonlinearities, extending the NLSE and the DNLSE demonstrates clearly that suitable combi- 
nations of nonlinear terms may be added to the original integrable equations without changing 
their linear dispersion relation, and at the same time preserving their integrability. 

Moreover, one can apply the whole procedure to each of the k = 3,4,... equations in the 
known integrable NLSE and DNLSE hierarchies , containing kth order linear dispersion. Each of 
such equations therefore would generate its new integrable EK hierarchy with higher nonlinear 
terms, but with the same linear dispersive term. 

It may be recalled that the generalized higher NLSE appearing in many physical problems 
[2, 3, 4], which contains higher nonlinear terms together with the same second-order linear 
dispersive term, can be reduced for particular parameter choices to integrable equations (7) or 
(8). Therefore we may hope that equations with hopelessly complicated nonlinear terms, arising 
in many other physical situations, could also be reducible to some of the higher nonlinear 
equations in the rich integrable EK hierarchies found here and hence could be exactly solved, at 
least for certain parameter choice or for some limiting values. The method presented here for 
constructing new integrable nonlinear hierarchies with unchanged linear dispersion is general 
enough to be applicable to other integrable PDEs with complex fields , e.g. mixed DNLS, 
complex mKdV etc. and even to vector models like Manakov model. 
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